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Current coronal mass ejection (CME) models set their lower boundary to be in the lower corona. 
They do not calculate accurately the transfer of free magnetic energy from the convection zone to 
the magnetically dominated corona because they model the effects of flux emergence using kine- 
matic boundary conditions or simply assume the appearance of flux at these heights. We test the 
importance of including dynamical flux emergence in CME modeling by simulating, in 2.5D, the 
emergence of sub-surface flux tubes into different coronal magnetic field configurations. We inves- 
tigate how much free magnetic energy, in the form of shear magnetic field, is transported from the 
convection zone to the corona, and whether dynamical flux emergence can drive CMEs. We find 
that multiple coronal flux ropes can be formed during flux emergence, and although they carry some 
shear field into the corona, the majority of shear field is confined to the lower atmosphere. Less 
than 10% of the magnetic energy in the corona is in the shear field, and this, combined with the fact 
that the coronal flux ropes bring up significant dense material, means that they do not erupt. Our 
r**"* ' results have significant implications for all CME models which rely on the transfer of free magnetic 

energy from the lower atmosphere into the corona but which do not explicitly model this transfer. 
Such studies of flux emergence and CMEs are timely, as we have new capabilities to observe this 
with Hinode and SDO, and therefore to test the models against observations. 



/\ ' Subject headings: CMEs, Flux Emergence, MHD 

.w. 

1. INTRODUCTION 

1.1. CME modeling 

Coronal mass ejections (CMEs) and eruptive flares are the most energetic manifestations of solar ac- 
tivity, with a typical CME accele rating 10 16 g of coronal plasma to speeds sometimes exceeding 1000 km/s 



([Gosling et al.lll976t lHildnerlll977l ). It is now generally accepted that these giant eruptions of solar plasma and 
field are due to the explosive release of magnetic energy stored in the corona prior to the event. Furthermore, 
the free energ y for CMEs/fl a res is believed to be stored in the strongly s heared flux of filament channels (see, 



e.g, reviews bv lForbesI J200QI) ; iKlimchukl (J200l( ); lLinton fc Moldwinl (J2009h ). This flux is sheared in that that the 



field lines are mainly parallel to the photospheric neutral line, rather than perpendicular as would be expected 
for a potential field. All CMEs/flares are associated with filament channels and these are the only locations 
in the corona where the magnetic field exhibits strong non-potentiality. The basic picture of a CME is that, 



prior to the event, there is a force balance in the corona between the upward magnetic pressure of the sheared 
filament channel field and the downward tension of overlying unsheared field. As the shear builds up or the 
overlying field evolves, this force balance eventually disrupts producing an explosive outward expansion of the 
filament channel flux and some of the overlying flux. 

It is evident from this discussion that the two central issues in understanding CMEs/cruptive flares are 
the process by which magnetic shear forms in filament channels, and the mechanism that disrupts the force 
balance. These two issues have been at the heart of CME and flare research for decades and are the features 
that distinguish the various CME models. The first issue of the shear formation process is especially critical, 
because the disruption mechanism is almost certain to depend on the topology of the filament channel field, 
which must be determined by the formation process. For example, some models invoke kink-like instabilities, 
which requires the formation of magnetic twist as well as shear in the filament channel field. We emphasize, 
however, that all the models and all the observations agree that the bulk of the magnetic free energy is in the 
form of shear. Note also that the issue of shear formation has broad implication for a number of important solar 
physics problems, in particular the origin of prominences and filaments. 

There are only two general processes that can produce the observed shear in filament channels: either the 
field emerges from below in a non-potential, sheared state, or the photospheric motions shear the field after it 
emerges into the corona. Of course, both processes must be present to some extent in the real corona. The 
latter process is likely to dominate in non-active region filament channels, such as those associated with high- 
latitude quiescent filaments, because these form well after the flux has emerged. On the other hand, in active 
regions, especially in the strong complex regions that are the source of fast CMEs, the filament channels form 
with the active region and, hence, the shear must emerge with the flux. Given that the shear formation is the 
fundamental driver of the eruption, it is clear that flux emergence must be explicitly included in any model for 
fast CMEs in order for the model to be physically rigorous. 

One such model, which generates fast CMEs from sheared filaments, but which should be modified to 



include the creation of thes e sheared filaments via flux emergence, is the magnetic breakout model (jAntiochos 



1998t lAntiochos et al.lll999l) . This model exploits the finite resistivity in the solar atmosphere to disrupt the pre- 



CME force balance by generating reconnection external to the filament channel. This reconnection reconfigures 
the magnetic field and allows the eruption of the highly sheared structure. The breakout model thus requires 
a multipolar flux distribution with a separatrix where reconnection can occur. Numerical simulations have 
investigated this model in both 2.5D, whe re the domain is axisymmetric and two dimensional but where al l 
three components of vector variables evolve ( Lynch et al.ll2 004; Ma cNeice et al.ll20 04; Dc Vore fc Antiochosll2005 ) 



and 3D, where the domain is three dimensional ( Lynch et al.ll2008tlDeVore fc Antioch os 2008). The basic model 



uses a current-free quadrupole magnetic field configuration, consisting of a central active region dipole embedded 
in an anti-parall el global dipole fiel d, with a null point separating the two systems. This is shown in Figure 



[TJa), taken from iLvnch et al.l (|2004f) . Figure HJa) shows, for this experiment, the vertical force balance in the 
central arcade along the latitude = 0° symmetry line for the initial equilibrium. The gravitational (pg) and 
gas pressure (— VP) forces balance, as do the magnetic tension (B.VB/47r) and magnetic pressure (-VB 2 /8ir). 
Shear flows are applied at the lower boundary, within the central arcade, to generate magnetic field parallel to 
the active region neutral line (in our 2.5D simulations the shear field will be the magnetic field in the ignorable 
direction). Figure [^b) shows the vertical forces after a short period of shearing. The presence of shear field 
exerts an upward magnetic pressure on the overlying field and stretches the original null point into a current 
sheet, as shown in Figures [TJb) and (c). Reconnection occurs when the current sheet collapses below grid scale 
and overlying field is removed by this reconnection, which allows further expansion of the central arcade. The 
outward expansion of the sheared arcade occurs exponentially as more and more overlying field is reconnected. 
Flare-like reconnection underneath the sheared dipole creates a disconnected flux rope which can then escape, 
as seen in Figures [ljd)-(f), at speeds comparable to observations of a fast CME. 

There are numerous other CME models which either assume the appearance of shea red magnetic flux at the 
low corona, or use drivin g velocities to create sheared structures from coronal equilibria. IChen fc Shibatal ( 20001 ) 



and lDubev et al.l ( 20061 ) simulated a 2.5D CME model where a pre-existing coronal flux rope, constrained by 



overlying field, is allowed to escape when kinematically imposed emerging flux causes reconnection in a filament 



channel beneath the flux rope. They also found that CME eruptions occur when kinem atic flux emergence 
is used to drive reconnection outside the filament channel. The flux cancellation model ( Amari et al.l 12000. 
20031 120071 ) is another model based on reconnection, but this model uses both boundary imposed shearing and 
the cancellation of magnetic flux at the neutral line of filament sites to initiate a CME. iFan fe Gibson! (J2007h 
developed a CME model where the lower corona is driven by the kinematic emergence of a sheared flux rope 
into a pre-existing coronal field. Here, loss of equilibrium via the kink or torus instability leads to an eruption, 
with the strength of the overlying field being an important constraint on the eruption. The breakout model has 
also been extended to cov er different magnetic field configurations, and also to include a solar wind, in both 
2.5D and 3D simulations |van der Hoist et al.l EoOtL IZuccarello et alJliool ISoenen et alj|2009h . Other CME 
models exist whic h do not directly rely on flux emergence, such as those based on the loss of equilibrium of a 
coronal flux rope ( Titov fc Demoulinlll999t iRoussev et al.ll2003 ; Kliem fc T orok 2006), and simulations of the 
rise of a pre-formed coronal flux rope into the solar wind by I Zhang fc Wul ( 2009 ). Although not relying on flux 
emergence directly, these models do assume that the flux ropes have already been formed in the corona and 
reached an equilibrium with the background coronal field, a state which is most likely caused by the emergence 
of magnetic field from beneath the surface. 

All these models, constrained by the need to extend the simulation domain to at least a few solar radii, 
do not model the lower solar atmosphere. The lower boundary of these simulations has a typical density of 
3 x 10 -16 g/cm and a typical plasma-/? < 0.2 where 



= 



£2/871"' 



(1) 



Based on the VAL model of the solar atmosphere (jVernazza et al.l Il98lh the photosphere, in contrast, has a 
density of around 3 x 10~ 7 g/cm and (3 > 1. These CME initiation models simply assume the appearance of 
flux at the low density, low j3 corona and do not self-consistently calculate a process for the flux emergence. 
These simulations do not, therefore, address the critical question of whether and how newly emerging, sheared 
magnetic flux can rise from its origins in the high (3 convection zone to the low corona where it is required to 
drive CME models. This is an important question to address now, as Hinode and SDO scientists are currently 
making detailed observations of flux emergence with vector magnetograms, which will allow the comparison of 
observations and state of the art CME models. 



1.2. Flux Emergence 



The current state of the theory of the formation of active regions suggests that dynamo actions in the solar 
interior create twisted, buoyant flux tubes which rise throug h the convec t ion zone, intersect the photosphere 



and create the obs e rved sunspots and bipolar active regions (|Parkerlll955t IZirirJll970t ISchuesslerlll979t IParker 



1988c IZwaanlll987t lEmonet fc Moreno-Insertislll998t iFisher et al.l 12000 ). The subsequent expansion of these 
flux tubes into the corona has been a subject of vigorous research for the past 30 years. The most likely 
driver of flux emergence into the lower corona is an ideal MHD instability driven by magnetic buoyancy (the 



magnetic buoyancy instability), a particular mode of which was first suggested bv IShibata et al 



2D simulations using a crude model atmosphere and magnetic field configuration. They found that the rise of 
magnetic flux into the corona due to this instability matched observed rise velocities of magnetic flux. 

The process of flux emergence has been studied more recently in 3D simulations with a more detailed 
model atmosphere that includes a convection zone, isothermal photosphere /chromosp here, transition region, and 
isoth e rmal corona, with a buoyant twist ed flux tube embedded in the convection zone ( Fanll2001tlManchester et al 



2004; 



Archontis et al.l |2004 iFanl 12009). The reader is directed to a comprehensive review of this subject by 
Archontisl ( 20081) . These simulations show that as the flux tube intersects the photosphere, bipolar active re- 



gions are formed which exhibit shear flows along their neutral line. As the tube continues to emerge into the 
corona, sigmoid-like structures are formed consistent with observations, with magnetic field being created par- 
allel to the neutral line (shear field). This suggests that flux emergence may be able to provide magnetic shear 
field sufficient to drive CMEs, but it is unclear how this magnetic shear emerges into the corona, as the axis of 



the emerging flux tubes remains rooted near the surface. 



Manchester et al.l (|2004l ) showed that shearing motions associated with the emergence are able to create 
a sheared coronal flux rope which rises out of the simulation domain, although the authors do not associate 
this with an eruption because the rise speed is limited by the amount of available 'axial flux' to drive an 
eruption. Note that in this paper we use the term flux tube to refer to the sub-surface tubes which we initiate 
our simulations with, and flux rope to refer to the twisted structures formed in the corona. 3D simulations 
(|Fanl 120091 ) have shown that without a pre-existing coronal field, the rise of coronal flux ropes created during 
flux emergence is constrain ed by the expanding emer ging field of the original flux tube , or envelope field. In a 



similar manner, the work of 



' Archontis fc Hoodl (|2008l ) and lMacTaggart fe Hoodl (|2009aD showed that flux ropes 
formed from multiple flux tube emergence are ultimately trapped by the envelope fields. For an eruption to 
occur, the overlying tension of the envelope field must be removed. 



Archontis fc Torokl (|2008[ ) added a pre-existing horizontal coronal field to their flux emergence simulations. 
They found that a coronal flux rope could be formed by inflows and horizontal shearing motions which drive 
reconnection within the expanding volume of the flux tube. They also found that reconnection between the 
envelope field and the pre-existing field r educed the overlying tension of the envelope field and allowed the 
flux rope to erupt. iMacTaggart fc Hoodl (J2009H ) found a similar result, but used a toroidal sub-surface flux 
tube which was already arched towards th e surface, which allowed the flux tube axis to emerge bodily into 



the atmosphere. lArchontis fc Hoodl (|2010l ) showed that the likelihood of a coronal flux rope escaping when 



reconnection occurs between the ambient coronal field and the envelope field depends strongly on the strength 
of the ambient field. 

Flux emergence simulations such as these occur over timescales of 1000 s, much quicker than timescales 
of CME models, and do not extend high enough in the corona to allow a full comparison with CME models 
and observations. The development of coupled models which span the timescales and spatial scales of both flux 
emergence and CME eruption is very important to fully understand the initiation of these explosive events. 

In this work we investigate whether magnetic flux emergence, in a 2.5D cartesian geometry, can provide 
enough magnetic shear energy to drive CME initiation, with focus on the initiation mechanism of the magnetic 
breakout model. Figure EJc) shows the energy in the shear magnetic field (normalized to the energy in the total 
field) in a typical breakout model simulation as shearing is applied at the surface. Also shown is the total kinetic 
energy in the system. As the inner arcade is sheared at the /3 = 0.2 lower boundary the shear energy increases. 
When the magnetic shear energy is about 20% of the total magnetic energy, the upward magnetic pressure due to 
the shear field can force reconnection at the null point and hence further expansion, resulting in a rapid increase 
in kinetic energy. Hence the breakout model typically needs 20% of the magnetic energy above the /3 = 0.2 
corona to be in the shear field. Although we know that flux emergence can create the shear flows and sheared 
magnetic fields near the photosphere which are similar to those required in the low corona by CME models, 
it has not yet been calculated how much shear energy the emergence of magnetic flux can actually supply to 
the low /3 corona. This issue is key to understanding how to couple dynamical flux emergence simulations and 
current CME simulations. We will investigate this using simulations of the evolution of sub-surface flux tubes 
with a variety of tube profiles as they emerge into both field-free corona and corona with pre-existing fields. We 
will also investigate whether flux emergence can drive CME initiation in these configurations. 

This paper is structured as follows: §2 explains the model, including numerical method, atmospheric and 
magnetic models, and initialization. §3 presents results from 2.5D simulations. Firstly we look at emergence 
into a field-free corona, then we look at emergence into a quadrupole, and finally emergence into a dipole. In 
all cases we attempt to cover a realistic range of field profiles and coronal field strengths to test the robustness 
of our results. §4 summarizes our findings and suggests consequences for modeling the initiation of CMEs by 
flux emergence. 



MODEL DESCRIPTION 



2.1. Numerical Method 



We model our system of a magnetized plasma with a polytropic equation of state (polytropic index 7 = 5/3) 
using the resistive MHD equations. These are solved numerically using the Lagrangian remap code LareXd 
(jArber et al.ll2001l ). The equations solved are presented below in Lagrangian form, using Gaussian units: 
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The gas density, pressure, and internal specific energy density are denoted by p, P, and e respectively and 
are defined at the center of each numerical cell. The magnetic field, denoted by B, is defined at cell faces, 
j = cV A B/4-7T is the current density, and c is the speed of light. The velocity, v, is defined at cell vertices. 
This staggered grid preserves V.B during the simulation. The gravitational acceleration is denoted by g, v is 
the viscosity, set to 2.6 x 10 4 g.cm _1 s _1 , r\ is the resistivity, and S is the stress tenor which has components 

representing an ideal gas: P — pRT, where R is the gas constant. 



i<5ijV.v), with Qj = ^(Iff 1 + tt^-)- The four equation are closed with a simple equation of state 



Equation [5j mi nus the last ter m, (e — eo(/?))/r, is similar to that used in the majority of flux emergence 
simulations to date ( Archon tis 2008). When flux tubes emerge and expand into the coron a, the associated coolin g 
due to the pressure ter ms can lead to unrealistic temperatures in the corona, as shown bv lLeake fc Arberl (|2006f ). 
Leake fc Arberl (J2006I ) also showed that by including a relaxation term, such as the last term in Equation [5j 
this unrealistic cooling could be avoided. This approach is a simple way to model the effects of terms that we 
are currently unable to formulate numerical equations for (such as chromospheric/coronal heating), and terms 
such as thermal conduction and radiative losses which are beyond the scope of this particular numerical code. 
The term (e — cq(p))/t relaxes the specific internal energy density back to its initial equilibrium values eo(p) on 
a timescale r which scales with density: 



r(y) = 23 



p(y) 
p(y = 0) 



-1.43 



(0) 



The equations are solved in 2.5D: the simulation box is 2D, with x and y being independent variables and 
z being ignorable, but all three components of the vector variables are evolved. The simulation box extends 
from -3 Mm to 90 Mm in the vertical direction (y) and from -45 Mm to 45 Mm in the horizontal d i rection (x) . 
This is much larger than most flux emergence simulations typically cover ( Manchester et al.ll2004t iFanl 120091) . 
but is designed to allow expansion of the flux tubes into coronal field as far as possible. The numerical grid is 
stretched in both x and y to provide better resolution in areas of interest. The grid spacing, S, is smallest at the 
surface, where S — 0.0375 Mm, 0.25 times the scale height at the surface, and largest at the upper boundary, 
where 5 = 0.15 Mm. 

In the majority of previous flux emergence simulations either ideal MHD is used {rj = 0), or 77 is chosen 
to be a constant such that the diffusion due to this explicit resistivity is greater than the numerical diffusion 
in the code. Theoretical evidence suggests that when the electron fluid flow speed, v e , exceeds the phase speed 
of the ion-acoustic mod el, c, , then ion-acoustic turbulence has a strong effect on current sheet development 
(|Bvchenkov et al.lll988l) . In this case an effective formula for anomalous resistivity would be 



■q = 7/0 max ( 0, I 

Cjn, 



(7) 



which we can rewrite, using v e ~ |j|/«e, as 



77 = 770 max 0, 



- 1 



Jcrit 



(8) 



where jcrit = neci a . We choose 770 = 1.6 x 10 _ s, and jcrit — 2.07 x 10 statampere.cm -2 , and this anomalous 
diffusion exceeds the numeric al diffusion in the code. This approach has been used in MHD simulations of kink 
instabilities of coronal loops ( Arber et al.lll999t iGerrard et al.ll2001 ). 



2.2. Initial Conditions 

2.2.1. Background Stratification 

The initial stratification is a simple ID model of the temperature profile of the Sun which includes the upper 
3 Mm of the convection zone, a photosphere/chromosphere, transition zone, and corona. The temperature profile 
is given by 

' T ph - T^y, y < 
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The convection zone profile (y < 0) is a linear polytrope which is marginally unstable to convection, with 
m = —zj being the adiabatic index for a polytrope. The temperature in the photosphere and chromosphere is 
assumed to be constant, T p f t — 5700 K, as is the temperature in the corona, T cor = 8.6 x 10 5 K. The height 
of the transition region is y tr = 3.75 Mm, and its width is w tr =0.75 Mm. The density and gas pressure are 
specified by initially assuming hydrostatic equilibrium. The resulting stratification is shown in Figure [21a) . 



2.2.2. Magnetic Flux Tube 

We insert a cylindrical magnetic flux tube into the model convection zone at a height of y t 
The axial field for all the tubes is given by 

B z = B e- r ' a 



-1.8 Mm. 



(10) 



where r = y/x 2 + (y — y t ) 2 is the radial distance from the center of the tube, B = 7800 G is the axial field 
strength at the center (r = 0), and a — 0.3 Mm is the width of the tube. This axial field will equate to the shear 
field required by the breakout model. We then choose three different tubes by choosing three different profiles 
for the twist magnetic field (Bg in the cylindrical coordinate system, B x and B y in cartesian). A minimum 
amount of twist is required to play the role of sur face tension and keep the flux tube coherent as it rises in the 
convection zone ( Emonet fc Moreno-Insertislll998 ). We denote the twist by 9, defined as 
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The three twist profiles we use are given by 
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where c = 1/a, c\ = 0.15 Mm and C2 = 1.5 Mm. 



Tube 1, the constant twist tube, is similar to tubes which have been studied exten sively in both 2D and 
3D (|Fanll200lirManchester et al.ll2004 Urber et al1l2007t iLeake fc Arberll2006l ; lFanll2009[) . This initial condition 
has been shown to create a buoyant flux rope that can not only emerge into the solar atmosphere, but interact 



with a simple, horizontal pre-existing coronal field and create isolated flux ropes in the corona, and is thus an 
ideal choice for our study. 



Tube 2 (increasing twist) and Tube 3 (decreasing twist), are similar to two tubes investigated bv lMurrav fc Hood 



(2008) and have twist profiles which increase and decrease with radius, respectively. In their simulations inves- 
tigating the effect of non-constant twist, the tops of these two flux tubes emerged higher into the atmosphere 
than tubes with other twist profiles and are therefore the most relevant to this study. 



As in iMurrav fc Hoodl ( 20081) . we can group the three tubes into low magnetic tension and high magnetic 
tension cases. The magnetic tension for these tubes is given by 



(B.V)B 
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(15) 



The magnetic field, twist (0), and tension profiles for the three tubes are shown in Figure|3J Tubes 1 and 3 have 
similar twist and tension profiles. Tube 2 has lower te nsion near its axis due to the lower magnitude of Bg close 
to the center of the tube. From their 3D simulations, IMurrav fc Hoodl (|2008l ) concluded that the twist profile 
was not as important a factor in the emergence of the flux tubes as the axial field strength, which affected the 
rate of emergence but not the extent of the emergence. The purpose of choosing these three tubes in our study 
is to cover a range of flux tube behavior whilst keeping the potential amount of shear energy provided by flux 
emergence optimal. 

The plasma (3 at the center of the flux tubes is approximately 10, and the flux tubes are all non-force- 
frcc. To initiate a buoyant flux tube we perturb the background gas pressure by an amount p\ such that 
(Vpi) = (j A B) , which can be integrated to give 
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B 2 {r) 
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dr 



(16) 



so that the tube is in radial force balance. Assuming that the flux tube is in thermal equilibrium with its 
surroundings makes the tube less dense that the surrounding plasma and initiates its buoyant rise to the 
surface. 

We will investigate the emergence of these tubes first into a field-free corona and then into different 
quadrupole and dipole coronal field configurations to test how efficient flux emergence is at providing mag- 
netic shear energy to the corona. 



2.2.3. Quadrupole Background Magnetic Field 

After investigating the emergence of flux tubes into a field free corona, we investigate the interaction 
of coronal field and emerging flux within the context of magnetic breakout. To do t his we first impose a 
quadrupole coronal field above the magnetic flux tube. We construct the quadrupole as in Karpen et al.l (J1996I ) 
by first defining the magnetic vector potential along a source surface at y = — 7.5Mm (which is outside the 
computational domain) : 



A x (x,-7.5 Mm) 




i - (x/x a y 



< \x\ < x a 

\x\ > X a 



(17) 



with x a = 21 Mm being the horizontal extent of vertical qaudrupole field at this height. The vector potential 
in the interior of the domain is then determined by integration of Laplace's equation with Equation \T7\ as the 
lower boundary condition and A z = at y = CO. The quadrupole consists of an inner dipole and an overlying 
dipole of opposite orientation, separated by a null point at y ~ 3 Mm. The height of the null is chosen to allow 
for the expansion of the flu x tube as it emerge s into the lower corona, during which the cross-section of the tube 
increases by a factor of 10 (jArber et al J 120071 ). Thus the 0.3 Mm flux tube will stay within the central arcade 
of the quadrupole when it emerges. The magnetic configuration of the sub-surface flux tube and overlying 



quadrupole field is shown in Figure [5] For the breakout mechanism, most of the horizontal quadrupolc flux 
above the null is removed by reconnection with horizontal flux below the null. We therefore vary the strength 
of the quadrupole so that the horizontal flux above the null point, along the x = line, is a certain fraction 
of the horizontal flux contained between the center of the sub-surface flux tube and the null point. We choose 
<7d to be [11.7, 5.85, 3.9] xlO 7 Gem, which give quadrupole surface strengths at X = x a /2 of [197, 98, 66] G 
respectively. The horizontal flux above the null point is then [1, 0.5, 0.3] times the horizontal flux between the 
flux tube and the null point. This allows us to control the amount of reconnection at the null point, whilst 
still being able to fix the at certain heights above the surface. The curves in Figure EJb) shows the profiles 
for the configurations of a sub-surface flux tube plus three different stre ngth quadrupoles. These profil es are 
genera lly consistent with profile models of active regions developed by iGarv fe Alexander! (|1999T) and I Gary 
POOII ) which take into account Soft X-ray Telescope limb observations. 



2.2.4-. Dipole BackGround Magnetic Field 

While emerging a flux tube into a quadrupole is the first step to coupling flux emergence and the breakout 
model, we also investigate a more self-consistent configuration, with a flux tube emerging into a background 
dipole field. By using a dipole with opposite orientation to the upper half of the flux tube, reconnection between 
these two systems creates a quadrupole structure self-consistently. 



The dipole field is represented by the vector potential A — A z e z where 



A z (x,y) = d 



y-Vd 



(18) 



with r\ = ^Jx 1 + (y — yd) 2 being the distance from the source. We chose yd to be -15 Mm so that the flux 
tube is far from the source of the dipole. To cover various dipole strengths we pick a range of d = [18, 9, 
4.5] xlO 28 Gem 3 which gives a magnetic field strength at x — 0,y = of [105, 53, 26] G respectively. In a similar 
fashion to the quadrupole field, we choose these values so that the horizontal flux contained between the axis 
of the flux tube and the separatrix which separates the flux tube and the dipole is a certain fraction of the 
horizontal flux above this separatrix. With decreasing dipole field strength, these factors are [1, 0.5, 0.25]. The 
profile for the three dipole choices overlying a sub-surface flux tube are shown in Figure [3jb). These choices 
of dipole strength allow f or a range in the profile s, bu t are s till co nsistent with the models of in the solar 
atmosphere developed bv lGarv fe Alexander! (J1999I ) and lGarvl (|200lh . 



RESULTS 



3.1. Initial Evolution in the Convection Zone and Lower Atmosphere: Effect of Twist Profiles 

The evolution of the flux tubes in the convection zone is similar to previous simulations of rising flux tubes 
( Emonet fc Moreno-Insertislll998l : lMagaral2001tlLeake fc Arberll2006l ). The initial buoyancy causes the flux tube 
to rise to the surface, during which time the cross-section increases and flux conservation decreases the axial, or 
shear, field strength. As the flux tube meets the convectively stable photosphere its rise is halted and horizontal 
expansion spreads the flux tube out to form a contact layer with the plasma above. The density above this 
layer is higher than the density below due to the concentration of magnetic field and the layer is unstable to a 
Rayleigh- Taylor like instability known as the magnetic buoyancy instability ( Newcomblll96lUAthav fc Thomas! 
1961uAchesonlll979 ). The stability of the contact layer involves a competitio n between the destabilizing gradient 
in the magnetic field and the stabilizing sub-adiabatic temperature gradient ( Gilmanlll970l h The stabilizing term 
is dependent on the local 0, and as the center of the tube rises up to this contact layer, the local falls and 
the in stability allows the upper portion of the flux tube to expand into the atmosphere, as in lArchontis et al.l 
d2004h . 



As can be seen in Figure Etc), the choice of twist profile affects the amount of tension in the tubes. The 



increasing twist tube (Tube 2) has relatively less tension near its axis than the constant twist tube (Tube 1) 
and the decreasing twist tube (Tube 3), and also a higher buoyancy at its center. This difference distinguishes 
the evolution of the increasing twist tube from the other two, and from now on we refer to the low tension 
case to mean Tube 2 and the high tension case to mean Tubes 1 and 3. In Figure [5] we can see how the choice 
of twist profile affects the initial expansion into the field-free corona. The tubes act nearly identically as they 
rise buoyantly in the convection zone and also in the lower atmosphere, as the magnetic buoyancy instability 
develops and the outer fieldlines extend into the corona. However, the two cases then differ in how the center 
of the tube reacts. In the high tension cases, the lower fieldlines near the center of the tube remain near the 
surface, shown in Figure [Hid), and the shear field is concentrated at the center, creating a single neutral line. In 
the low tension case, the decreased tension allows the center of the flux tube to rise, forming a crescent shape, 
shown in Figure EJh). During this evolution, the mass in the center of the tube drains to the lowest possible 
location which is at either end of this crescent shape. The coupling of shear field and density in this 2.5D 
simulation means that the shear field also concentrates in these regions. Note that this emerging structure now 
has three neutral lines at the photosphere. 



3.2. Emergence into a Field-Free Corona: Effect of Twist Profiles 

The center of the high tension flux tube remains at the surface for the duration of the simulation. We 
follow the continued evolution of the low tension flux tube in Figure [7J which shows the shear field and the 
log of density, along with magnetic fieldlines. The deformation of the center of the tube into a crescent shape 
forms a current sheet as field of opposite sign is brought together on either side of this crescent. This enhanced 
current triggers the anomalous resistivity and reconnection sets in. As can be seen from Figures [7jb) and (e), 
the reconnection sites are at either side of the apex of the current sheet arch. Plasmoids are formed which 
are expelled along the current sheet. The half width of the current sheet, I, is approximately 0.15 Mm, and 
a typical distance between the plasmoids is 5 Mm. Thu s the ratio of ha l f curr ent sheet width to the length 



between islands, A, is approximately l/X — 0.03. Based on lPriest fc Forbes! (|2000h . the wavelength range for the 
tearing mode is 

l(D\ 1/A I 



2,\lv A J < X < 1 /^)-0.16 (19) 

where D is the magnetic diffusivity (D — r/c 2 /Ait) and va is the Alfven speed. Using values from our simulations 
gives the lower limit as 0.01. Hence the current sheet dimensions, l/X — 0.03, are compatible with a tearing 
mode instability. Later in time, plasmoids are formed closer to the apex of the current sheet and coalesce, 
creating an isolated flux rope at the center. Both shear field and mass are trapped in this flux rope, but the 
majority of the shear field concentrates at either end of the current sheet, in the low-lying dips in the magnetic 
field. 

To determine what happens to this coronal flux rope which is formed from the original center of the low 
tension flux tube, we can examine the forces in the flux rope. Figure [8] shows the four vertical forces associated 
with gravity, magnetic tension, magnetic pressure, and gas pressure along the central, x = 0, line at two 
different times. The forces are shown for the initial sub-surface flux tubes in Figures Eta) and (b), where the 
initial hydrostatic equilibrium can be seen in the large scale trends for pg and — VP, and the magnetic forces 
can be seen over the range -2.5 Mm to -1.5 Mm. The outward directed magnetic pressure force, due primarily 
to the shear field, dominates over the inward directed tension force of the twist field. The gas pressure force 
is directed inwards to balance this excess magnetic pressure force. Figures [S[c) and (d) show the forces in the 
corona after the flux tube has emerged. For the low tension case, the forces in the newly formed coronal flux 
rope, the center of which is at y = 7 Mm, can be seen in Figure [H^d). The presence of mass can be detected in 
the small dip in the gravitational force. Here, in contrast to the configuration of the initial flux tube, the inward 
directed tension of the twist field dominates over the outward directed magnetic pressure force. Thus the gas 
pressure force is directed outward here, to balance the excess tension force. This shows that, in contrast to the 
breakout model, shear field does not create a strong unbalanced upward/outward pressure force, and so there is 
no impetus which can drive a CME eruption. This shows the dramatic effect the dense lower solar atmosphere 
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can have on the role of shear in driving CMEs. 

As far as we are aware, the mechanism shown here for the creation of a coronal flux rope from an emerging 
flux tube whose cross section distorts into a crescent shape has not been seen before. We expect that such 
coronal flux ropes will form from 3D flux tubes which emerge with a large radius of curvature along their axis. 
A large radius of curvature implies that the draining of mass along the axis will not exceed the draining of 
mass perpendicular to the axis seen in these 2.5D simulations, and so the behavior will be similar. However, 
3D simulations of similar flux tubes by iMurrav fc Hoodl ( 20081 ) do not show this behavior. We find that we 



still see the formation of our coronal flux ropes when we use the same parameters (yt =-0.17 Mm, J3o= 3 900 G , 



<z=0.425 Mm, c = 3 Mm and ci=0.17 Mm) as one of the low tension tubes used in lMurrav fc Hoodl (|2008i ). 
3D simulations are required to investigate this issue further. 

Figure [9ja) shows the height of the centers of the three flux tubes as they emerge into the field-free corona. 
The center (which is initially the axis) is defined at the point where B x \ x= q changes sign. The centers of the 
high tension tubes stay at the photosphere, while the center of the low tension tube rises to a height of about 
8Mm as the tube is deformed into a crescent shaped current sheet and the resulting coronal flux rope rises. The 
coronal flux rope does not escape, most likely due to the presence of mass accumulation in the rope and due to 
the lack of magnetic shear. 

Figures [9fb) and (c) show the normalized magnetic shear energy above two different heights as a function 
of time for all three tubes. This normalized energy is given by 



ry=90Mm f x=45Afm r2 , , 
Jy= yi Jx=-4:5Mr,iz uxa y 

,,,, \B\ z dxay 



77. _ J V=Vl Jx=-4bMm z o /r,„N 

Shear — rV=90Mm r x=45Mm i n i,j_j_ v zu J 



The first height, yi = 1.2 Mm, is the lowest height at which f3 = 0.2 in the simulations in this paper with initial 
non-zero coronal field, and is an estimate of the height at which we need to transfer a significant amount of the 
shear magnetic energy relative to the total magnetic energy in order to drive a CME. The second height is at 
y\ = 5 Mm above the surface and is at a location just below where the isolated flux rope is formed from the 
center of the low tension tube. 

As the low tension case develops (£=1700 s to £=2200 s), the crescent shaped current sheet is formed and 
mass and shear field drain to the ends, below 1.2 Mm. Thus the shear energy above 1.2 Mm drops as shown by 
the green/dashed curve in Figure [D^b) . The coalescence of plasmoids, which leads to a concentration of shear 
flux, then creates an increase in E s h ea r at £=2300 s. The expulsion of these plasmoids along the current sheet, 
down below 1.2 Mm, creates the sharp decrease seen shortly after at £=2400 s. As the flux rope formed from 
this current sheet then rises into the corona, carrying shear field with it, the shear energy above 5 Mm rises, as 
can be seen in the green/dashed plot in Figure EJc). The gradual rise in E s h ear for the high tension cases above 
1.2 Mm after £=2800 s (red/dot-dashed and blue/solid plots in Figure^b)) is caused by reconnection at nearly 
vertical current sheets which form at the edges of the active region. This reconnection creates plasmoids which 
are able to transfer a small amount of shear above 1.2 Mm but which do not rise to the 5 Mm level. 

The magnetic shear energy above 1.2 Mm for all three tubes does not exceed 10% of the total magnetic 
energy above 1.2 Mm. For the emergence of a 2.5D flux tube into a field-free corona, the shear field is generally 
confined to the high /3 lower atmosphere. The shear field which does rise in the coronal flux rope is not strong 
enough to create magnetic pressure which will counterbalance the magnetic tension and gravitational force in the 
rope and hence drive a further rise or an eruption. In the next section we look at the expansion into a quadrupole 
field, to see if the small amount of shear field which does emerge into the corona can drive reconnection at the 
null point of the quadrupole, thus reconfiguring the field to allow an eruption as in the breakout model. 



3.3. Emergence into a Quadrupole Coronal Field 

We present results of the emergence of the low tension flux tube into a pre-existing coronal quadrupole 
magnetic field. The quadrupole configuration can be destabilized by the presence of magnetic shear, as in the 
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breakout model, but the effects of the lower atmosphere have not yet been investigated. We know from the 
previous section that the emergence of magnetic flux can bring shear into the corona, but for a field-free corona, 
the coronal flux ropes which then form are not sheared enough to erupt. In this section we emerge the low 
tension flux tube into a quadrupole of varying strength, and examine the interaction between emerging field and 
coronal field. We also look at the effect of the quadrupole strength on the overall emergence. We test whether 
the quadrupole configuration can be destabilized by flux emergence-supplied magnetic shear energy. The initial 
configuration is shown in Figure [5]for the medium strength quadrupole (B\ x=Xa /2,y=o — 98 G). 

The initial rise of the flux tube is similar to the field-free corona case up to about 1000 s. The presence of 
quadrupole field below the surface has a negligible effect on the rise of the tube, and its initial expansion into the 
atmosphere. Figure [10] shows the evolution of the low tension flux tube and the medium strength quadrupole 
at 6 different times (the coupling of B z and p in these 2.5D simulations allows us to only show the density to 
demonstrate the structure of the coronal flux ropes) . The field of the inner dipole of the quadrupole is aligned 
with the twist field of the upper part of the tube which expands into the corona. Thus the inner flux system 
of the quadrupole field is pushed upwards towards the null by the emergence, as shown in Figures [TOfa) and 
(b) . The expanding shell of the tube causes a crescent shaped current sheet to be formed as the original null is 
stretched over this shell. Reconnection sets in on either side of the apex of the current sheet. The current sheet 
continues to rise as the expanding shell of the emerging flux tube pushes it up as shown in Figures [TUTd )- (f ) . 
while the reconnection creates a coronal flux rope at a height of 18 Mm. 

In Figure [10] we can also see the deformation of the emerging flux tube's center, and the creation of a flux 
rope by reconnection at the current sheet formed at about 5 Mm, just as in the field-free corona simulations. 
Thus there are two coronal flux ropes formed, as can be seen in Figure [TOTf ) . The upper rope is created by 
the reconnection at the separatrix between the two systems of the quadrupole, driven by the expansion of the 
flux tube into the atmosphere. The lower rope is formed by reconnection at the current sheet formed from the 
center of the emerging flux tube. Both of these flux ropes carry shear field and mass. 

Figure [IT] shows the vertical forces in these two coronal flux ropes at later times of t = 2392 s and 
t = 2760 s. Both ropes carry significant amounts of mass and in both ropes the magnetic tension acts inwards, 
counteracting any outward force due to gas or magnetic pressure. Figure lllf a) shows that both the magnetic 
tension and magnetic pressure forces are directed inwards in the upper flux rope, thus a strong outward pressure 
force is required for force balance. In addition, the downward gravity force of the trapped mass is as strong 
as any of the pressure or tension forces. These other forces all have an upward directed asymmetry to help 
balance this large gravity force, but nonetheless, as can be seen in Figure [UJc), they cannot and the flux rope 
starts to sink. In contrast, the gravity force is relatively small in the lower flux rope, as shown in lllf b). and the 
magnetic pressure force is directed outwards. Yet this outward magnetic pressure force is well balanced by the 
inward magnetic tension. Thus, while the ropes differ in the way the magnetic and gas pressures act, the key 
fact is that the amount of shear field in them is not enough to create an unbalanced outward acting magnetic 
pressure force to cause any further expansion. 

The fraction of the magnetic energy contained in the shear field above the heights of 1.2 Mm and 5 Mm for 
the three choices of quadrupole and the field- free corona case are shown in Figure [12] In general, increasing the 
strength of the quadrupole decreases this fraction as the total magnetic energy in the system increases. However 
this also increases the amount of reconnection at the separatrix of the quadrupole as the configuration contains 
more flux above and below the null. 

We have shown that the emergence of a flux tube into a quadrupole field does cause reconnection at the 
separatrix between the two flux systems in the quadrupole. However, the nature of the reconnection means 
that it does not remove significant amounts of horizontal field from above and below the original null point, 
as in the single X-point reconnection seen in the breakout model. The reconnection at the null is driven here 
by expansion of the outer part of the emerging flux tube, and the current sheet formed prefers reconnection 
off-center, creating a coronal flux rope which does not have enough magnetic pressure to escape. This flux rope 
eventually falls back down, and joins with the flux rope which was created from the center of the emerging low 
tension flux tube. The magnetic energy in the shear field above the j3 = 0.2 level is too small (less than 10% of 
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the total magnetic energy) to provide the shear energy required by the breakout model, which is consistent with 
the fact that the shear field in the flux ropes is not strong enough to overcome the tension and gravitational 
forces. 

In the next section we show the emergence of flux tubes into a dipole coronal field which is orientated 
opposite to the upper part of the emerging flux tube. Reconnection between the expanding tube and the pre- 
existing dipole coronal field may remove twist field in the coronal flux ropes formed and allow them to rise 
further. 



3.4. Emergence into a Dipole Coronal Field: Effect of Tension Profile 

In this section we present results from simulations of the emergence of the three different flux tubes into a 
particular dipole field to highlight the importance of the choice of tube profile. We choose the dipole field with 
•B| x =o,y=o = 52 G. This choice of dipole best highlights the differences between the high and low tension cases. 
We will investigate the choice of dipole strength in <j3.5l 



Figure [13] shows the emergence of the high and low tension tubes into the dipole field. The behavior of the 
flux tubes prior to emergence is very similar to the field-free corona and quadrupole coronal field cases. The 
outer part of the flux tube and the coronal dipole field lie on either side of a separatrix. As the flux tubes rise 
into the dipole field this separatrix is pushed upwards and becomes a current sheet, as shown in Figures [TBT a) 
and (d). In both cases the current sheet undergoes reconnection, consistent with the tearing mode instability, 
creating plasmoids on either side of the apex which are expelled along the current sheet. Eventually a coronal 
flux rope is formed by coalescence of plasmoids near the apex of the current sheet. This flux rope has an initial 
rise speed of 25 km/s. In the low tension case, the flux rope is formed more rapidly and at a higher point in the 
corona. Also, a flux rope is formed in the lower corona by the deformation of the center of the emerging low 
tension flux tube, leaving two flux ropes in the corona. The structure of these two flux ropes is similar to that 
of the flux ropes formed in the field-free and quadrupole coronal field experiment, with magnetic tension and 
a gravitational force acting to prevent any further rise or expansion. The creation of these coronal flux ropes 
appears to be the tearing mode instability, as w as the case in simulations of flux emergence into a pre-existing 



field performed by lArchontis et al.l ( 2006 L 2007 ). However this is a different mechanism from that which creates 



the flux ropes in the simulations of Manchester et al.l ( 2004 ) . lArchontis fc Torokl ( 20081) . and iFanl ( 2009 ). which 



is based on shear flows, converging motions, and internal reconnection of an emerging flux tube. 

Figure fMt a) shows the height of the flux tube centers (solid lines) and the height of the scparatrices between 
flux tube and dipole (dashed lines), along the x = line for the three tubes. Initially the paths of the centers 
and separatrices are the same for all three tubes as they expand and push the dipole field upwards. After 
t = 1500 s the low tension tube undergoes reconnection at its center, forming the lower of the two flux ropes. 
Reconnection at the separatrix between the emerging tube and the dipole field in all three simulations creates 
the coronal rope higher in the corona (solid lines). These coronal flux ropes eventually fall back down through 
the transition region after t = 4000 s. 

The amount of normalized magnetic shear energy supplied to the two heights of 1.2 Mm and 5 Mm for 
all three tubes emerging into the dipole can be seen in Figures RUT b) and Etc) . The normalized shear energy 
above 1.2 Mm increases after 2000 s for the experiment where the low tension tube emerges into the dipole, 
shown as the green/dashed curves on Figure [J"4Tb). This is caused by the removal of some horizontal flux by 
reconnection in the current sheets. While this is similar to the reconnection which leads to an eruption in the 
breakout model, this reconnection does not remove enough horizontal field to raise the importance of the shear 
field in the coronal flux ropes formed by this reconnection. As a result the shear magnetic energy is still less 
than 10% of the total magnetic energy above the (3 = 0.2 corona. 3D simulations bv lArchontis fc Torokl ( 20081 ) 



showed that reconnection between emerging flux tubes and horizontal coronal field allows coronal flux ropes 
which are formed during flux emergence to escape. This is not the case in these 2.5D simulations with a dipole 
coronal field. In the next section we investigate the importance of the dipole strength on the behavior of the 
coronal flux ropes. 
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3.5. Emergence into a Dipole Coronal Field: Effect of dipole field 

By varying the dipole strength we can control the location and amount of reconnection at the current 
sheet formed between the flux tube and the dipole, and through this affect the amount of shear energy (relative 
to total magnetic energy) which is transferred to the low /3 corona. In this section we present results of the 
emergence of the low tension flux tube into a range of dipoles. These dipoles have horizontal fluxes which are 
factors of the horizontal flux that is initially in the top half of the emerging tube of 1, 0.5 and 0.25. We also 
compare these three cases to the field-free corona case. 

The strength of the dipole affects the expansion of the outer part of the flux tube, as the stronger the dipole 
is the more tension it has to resist the expansion of the upper part of the flux tube into the corona. The stronger 
the dipole, the quicker (and lower down in the corona) the two systems are forced together by the emergence 
and so the quicker the reconnection occurs between them and the quicker the coronal flux rope is formed. This 
can be seen in Figure [T5] which shows snapshots at t=2530 s for the low tension flux tube emerging into the 
3 different dipoles (of decreasing strength from left to right) and a field- free corona (far right panel). At this 
time, the stronger dipole case has already formed a coronal flux rope, while the weaker dipole case has only just 
started reconnection between the emerging flux tube and the dipole. 

Figure ITST a) shows the heights of the flux tube centers and the flux tube/dipole separatrices for the different 
choices in dipole strength. The separatrices all reach similar heights, but do so at different times, due to the 
differing tension in the dipoles which resists the pushing up of the dipole by the expanding flux tube. 

Figures [ToTb) and [Trite) show the normalized shear energy above 1.2 Mm and 5 Mm. In general, the 
magnetic shear energy due to flux emergence is confined too low in the atmosphere to have a strong enough 
effect on the magnetic field in the corona. The shear energy is still always below 10% of the total energy above 
the low fi corona. Any shear field that is transported to the corona by coronal flux ropes is too small to affect 
the twist field at these heights. The coronal flux ropes formed are not sheared enough and contain significant 
amounts of dense material. In order to drive a CME by flux emergence in these simulations, we need a method 
for transferring more shear field to coronal heights while transferring less mass. 



4. CONCLUSIONS 

We have investigated the initiation of CMEs by dynamical flux emergence using 2.5D cartesian simulations 
of the emergence of a range of sub-surface flux tubes into a range of coronal magnetic field configurations. Most 
current CME models do not include the lower dense atmosphere in their models, but the boundary conditions 
for the shear energy and dynamical flows they use are most likely results of flux emergence through this region. 
In a first attempt to couple the modeling of CME initiation and flux emergence, we have taken coronal field 
configurations which are based on the magnetic breakout model. However, our results are generally applicable 
to any CME model which assumes the appearance of highly sheared flux at the low j3 corona. 

The evolution of fl ux tubes into a field-free corona can be separated into low and high tension cases, as in 
Murray fc Hoodl ( 20081 ). based on the relative amount of tension near the center of the flux tube. During the 



emergence of a low tension tube, an isolated flux rope was created by the deformation of the flux tube center into 
a current sheet and the reconnection of fieldlines on either side of the apex of this current sheet. Since twisted 



flux ropes are essential fo r several CME models (JTitov fc Demoulinlll999t iRoussev et al.ll2003t iKliem fc Torok 



20061; iFan fc Gibsonll20071 ) further investigations are warranted, especially in 3D, to take into account the effect 
of curvature along the tube's axis. The comparison of 3D simulations with observations of neutral lines during 
flux emergence will be able to verify if this behavior is occurring on the Sun. 

In an attempt to test whether the shear field brought into the corona by flux emergence could destabilize 
the quadrupole configuration of the magnetic breakout model, we added a pre-existing quadrupole coronal field 
above the emerging flux tubes. The expanding outer shell of the flux tube stretched the original null point of the 
quadrupole into a current sheet, and reconnection in this current sheet created another coronal flux rope, higher 
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in the corona. Both types of flux ropes were created in similar fashions: by the formation of a crescent shaped 
current sheet, with preferential reconnection occurring off-center, forming a coronal flux rope at the apex of the 
crescent. However, one type of flux rope was created by the deformation of the center of an emerging flux tube, 
and one was created along the separatrix between the two flux systems of the quadrupole. 

The emergence of sub-surface flux tubes into a simpler dipole coronal field configuration was simulated, 
in an attempt to remove the twist field of the flux tube which first emerges into the corona. Reconnection 
at the separatrix between emerging field and pre-existing dipole field yielded an isolated flux rope, as did the 
deformation of the center of the emerging low tension flux tube 

Although the coronal flux ropes that are formed in these simulations of flux emergence transport shear 
field into the low j3 corona, the shear field is weak compared to the twist field and they are laden with photo- 
spheric/chromospheric material. Any outward directed forces from the shear field are insufficient to overcome 
the inward directed tension force and the downward directed gravity force. These coronal flux ropes are therefore 
not candidates for erupting structures, and eventually disappear or fall back into the lower atmosphere. 

We have covered a range of flux tube profiles and coronal field structures and strengths, and in all cases 
found that the emergence of sub-surface flux tubes in these 2.5D simulations is an inefficient method for the 
transfer of magnetic shear energy into the low j3 corona. The amount of magnetic shear energy supplied to the 
low j3 corona was calculated, at heights of 1.2 Mm and 5 Mm. The lower of these two heights corresponds to the 
lowest point at which j3 = 0.2 in the simulations in this paper, and thus the lowest point at which the conditions 
are the same as those used in CME models. In all cases, the amount of magnetic shear energy supplied to the 
/3 < 0.2 corona was less than 10% of the total magnetic energy and appears insufficient to directly drive a CME. 
We conclude that the simple observation of flux emergence and shear at the photosphere is not sufficient evidence 
that a CME can be driven. Further investigations must be carried out to help understand the conditions under 
which emerging shear can be transported up to the low corona, where it is needed. 

For 3D flux tubes with significant curvature along the tube axis, flows along fieldlines may be an important 
phenomena in flux emergence simulations. These simulations should be repeated in 3D to see if this has an 
affect on the formation of coronal flux ropes and their eventual evolution in the solar atmosphere. 

We have tested the robustness of these results by varying the resolution and introducing asymmetry into 
the flux emergence and coronal fields, and find similar results regardless of the position of the flux tube relative 
to the center of the quadrupole/dipolc. 

In these simulations, and almost all previous flux emergence simulations, the high /3 lower atmosphere 
is assumed to be fully ionized, which is not the case for the solar chromosphere. The presence of neutrals 
in the nearly neutral regio n near the t e mper ature minimum may have an important effect on these results. 
Leake fc Arberl (2000) and lArber et al.l ( 20071 ) showed that this nearly neutral region can dramatically affect 



flux emergence. In particular, cross-field currents in emerging structures are destroyed by ion-neutral collisions, 
and the magnetic field lines are able to slip through the neutral gas during emergence. Thus the flux tubes lift 
up less mass from the chromosphere when the effects of partial ionization are taken into account. All this occurs 
lower down in the solar atmosphere around the /3 = 1 region and will directly affect the transfer of both mass 
and shear energy to the low ft corona above. These effects, along with 3D effects, may provide a mechanism for 
the transfer of magnetic shear energy into the corona without the transfer of mass, which may be key to driving 
a CME with dynamical flux emergence. This topic will be investigated in future publications. 
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Fig. 1. — Meridional projections of magnetic field l ines throughout the magnetic breakout eruption process at 
six different times. Reprinted with permission from lLvnch et al.l (2004J). Spatial units are solar radii. 



19 



E 1 5^ 
< (a) 



c 1.0 

T3 



0.5 



CO 
0) 

o 
o 



-pg 

+B'VB/4tt 

-V(B 2 )/8tt 

-VP 



a) 0.0 _-_-_-^-"-' 

u ; 



time = 0.00000 s 



E 1 5 

\ 

Q) 

c 1.0 

■D 



0.5 



o 
o 



1.05 1.10 1.15 1.20 
Solar Radii 



i 



-pg 

+B'VB/4tt 

\ -V(B 2 )/8tt 

\ -VP 



0.0 =~=r- -S ^t^~-: 



time.= 50000,0 s 



1.05 1.10 1.15 1.20 
Solar Radii 



0.6 pcT 
0.5 ! 

0.4 r 

0.3 \ 

0.2 r E 

0.1 t- 
0.0 L 




sheor ~//B 2 drd0 



40 60 

Time (1000 s) 



Fig. 2. — Panels (a) and (b): Vertical forces along the latitude=0° line at the center of the inner arcade of the 
magnetic quadrupole in the breakout model at t=0 s, panel (a), and t=30000 s, panel (b), after shearing of 
the arcade by imposed boundary flows. Panel (c): The normalized shear magnetic energy (solid line) and 2D 
integral of kinetic energy density (dashed line) in the system during this shearing period. 
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Fig. 3. — Panel (a): the initial height profile of the model stratified atmosphere showing the temperature (solid 
line), gas pressure (dashed line), and magnetic pressure (dot-dashed line), up to 20 Mm. The magnetic pressure 
profile is for a sub-surface magnetic flux tube and an overlying, anti-parallel dipole. Panel (b): The initial 
plasma ft for six different simulations: Three simulations with background dipole field and three simulations 
with background quadrupolc field. The straight black dotted line shows a constant value of ft = 0.2. 
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Fig. 4. — Panel (a): Initial twist function, O, for the three flux tubes. Tube 1: Constant twist, high tension, 
blue/solid lines. Tube 2: Increasing twist, low tension, green/dashed lines. Tube 3: Decreasing twist, high 
tension, red/dot-dashed lines. Panel (b): The axial field, B z , and twist field, Be, for the same three tubes. 
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Fig. 5. — Log of density at t = and initial magnetic field lines in the x — y plane (constant contour intervals 
of the z component of the magnetic vector potential, A z ). The field lines show the outer part of the sub-surface 
flux tube (located at x=0, y=-2 Mm), and an overlying quadrupole field. The inner arcade of the quadrupole 
is parallel to the upper half of the flux tube, as can be seen by the arrows, which show the direction of B x in 
the plane. The null is located at about 3.2 Mm above the surface (y=0). 
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Fig. 6. — Panels (a) to (d): The magnetic shear field, Bz (G), and magnetic field lines (constant contour intervals 
of A z ) at four different times for the high tension flux tube (Tube 1) emerging into a field- free corona. Panels 
(e) to (h): The same but for the low tension tube (Tube 2). 
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Fig. 7. — Panels (a) to (c): The magnetic shear field, Bz (G), and magnetic field lines, at three different times 
for the low tension flux tube emerging into a field-free corona. Panels (d) to (f): The same but for the log of 
the density. The two white tick marks on the left are located at y=l .2 Mm and 5 Mm. 
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Fig. 8. — Panels (a) and (b): Vertical forces along the central, x = 0, line inside the initial buoyant sub-surface 
flux tubes at t = s for the high tension and low tension cases. Panels (c) and (d): Vertical forces along the 
central, x = 0, line above the surface at t = 2760 s for the high tension and low tension cases. Panel (d) shows 
the coronal flux rope formed from the emerging low tension tube. 
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Fig. 9. — Panel (a): The height of the center of the three flux tubes emerging into a field-free corona. Panels (b) 
and (c): The normalized magnetic shear energy, -E^hear ( see Equationl20|). above 1.2 Mm and 5 Mm respectively, 
for the same three tubes. 
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Fig. 10. — Log of density and fieldlines for the low tension tube at six different times as it emerges into a 
medium strength (B\ x=Xa / 2 . y =o — 98 G) quadrupole field. The arrows show the direction of B x in the x — y 
plane. 
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Fig. 11. — Vertical forces along the central, x = 0, line for the two flux ropes formed when the low tension 
flux tube emerges into the medium strength (B| a . =Xa /2,i/=o = 98 G) quadrupole. Panel (a) shows the forces in 
the upper flux rope formed by reconnection at the current sheet formed from the initial quadrupole null point. 
Panel (b) shows the forces for the flux rope created by the deformation of the center of the emerging flux tube. 
Panels (c) and (d) show the total vertical forces in these two flux ropes, respectvely. 
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Fig. 12. — Normalized magnetic shear energy, E gnear (see Equation [20)) . above two different heights for the low 
tension flux tube emerging into a quadrupolar coronal field. Panel (a) is for 1.2 Mm. Panel (b) is for 5 Mm. 
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Fig. 13. — Log of density and fieldlines at three different times for the emergence of a high tension flux tube, 
panels (a) to (c), and a low tension flux tube, panels (d) to (f), as they emerge into the medium strength, 
B|x=o,y=o = 52 G, dipole coronal field. The arrows show the direction of the magnetic field. 
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Fig. 14. — Panel (a): Height of the flux tube centers (dotted lines) and tube/dipole separatrices (solid lines) 
for the three flux tubes emerging into a dipole of medium strength {B\ x=0 ^ y=0 — 52 G). Panels(b) and (c): The 
normalized magnetic shear energy, E g f lear (see Equation |20|) . above heights 1.2 Mm and 5 Mm for the same 
three experiments. 
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Fig. 15. — Shear field, B Z (G), and fieldlines at t=2530 s for the emergence of the low tension flux tube into 
4 coronal fields of decreasing strength. Panel a): B|x=o,2/=o = 104 G. Panel b): B\ x =o,y=0 — 52 G. Panel c): 
B|x=o,y=o — 26 G. Panel d): B\ x =o,y=o = G (no dipole). Panels (e) to (h) show the same but for the log of 
density. Arrows show direction of magnetic field. 
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Fig. 16. — Panel (a): Height of the flux tube centers, (dashed lines) and tube/dipole separatrices (solid lines) 
for the low tension flux tube emerging into dipoles of varying strength. Panels (b) and (c) : The normalized 
magnetic shear energy, £t near (see Equation 120)) . above heights 1.2 Mm and 5 Mm for the same experiments. 



